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Abstract. We report the results of our optical speckle interferometric observations of the nearby triple system 
GJ 795 performed with the 6 m BTA telescope with diffraction-limited angular resolution. The three components 
of the system were optically resolved for the first time. Position measurements allowed us to determine the elements 
of the inner orbit of the triple system. We use the measured magnitude differences to estimate the absolute mag- 
nitudes and spectral types of the components of the triple: M^ a =7.31±0.08, Af^ b =8.66±0.10, M^=8.42±0.10, 
SpAa «K5, SpAb ~K9, Spb «K8. The total mass of the system is equal to E_Mab=1.69±0.27.M©. We show 
GJ 795 to be a hierarchical triple system which satisfies the empirical stability criteria. 



1. INTRODUCTION 

According to the most recent concepts, stars form in small 
groups and clusters. The disruption of such groups results 
in the formation of both multiple systems and single stars. 
Stars in triple and more complex multiple systems make 
up for more than 20% of the Milky- Way population. The 
study of their dynamical and physical parameters is nec- 
essary for understanding the process of star formation as 
a whole. However, the now available observational data 
are insufficient for testing the theories of formation and 
evolution of multiple stars. We do not yet entirely under- 
stand the initial conditions and mechanisms of multiple 
star formation. The issues that still remain unclear include 
the conservation of angular momentum in the process of 
star formation; the dynamical stability of multiple systems 
with more than three stars; the effect of tidal interactions 
on the dynamical evolution of multiple systems; the distri- 
bution of orbital periods, eccentricities, component mass 
ratios, and correlations between these parameters; mutual 
orientations of the orbital planes in multiple systems, etc. 

Of special interest is the study of multiple systems 
with low degree of hierarchy with comparable orbital 
periods and semimajor axes. The Orion Trapezium — a 
small cluster of very young and massive stars — is the 
most well known dynamically unstable multiple system. 
The disruption time scale of the Trapezium is estimated 
at 10 4 - 10 6 years flPflamm-Altenburg fe Kroupa, 20061 ) . 
Main-sequence stars are also found in a number of systems 



that are potentially dynamically unstable QFekel, 1 981; 
|Szebehely fc Zare, 19771 |Tokovinin et al., 2005[ ). The 
stability criteria for multiple stellar systems were an- 
alyzed by Eggleton fc Kiseleva, 1995[ |Golubev, 1967 



Gol ubev, 19681 |Harrington,T972l |Harrington,T977] 
Steves & Roy, 1999 and other authors. However, the 



Send offprint requests to: E. V. Malogolovets, e-mail: 



conclusion about the dynamical instability of a particular 
system is often disproved when its orbital elements are 
refined. The triple system ADS 16904 with the periods 
of the inner and outer orbits equal to 15 and 150 
years, respectively, is an example (Bale ga et al., 19 99). 
According to all known criteria, the state of this system 
must be close to instability (Orlov & Zhuchkov, 2005). 
However, new interferometric measurements including 
the data obtained with the 6 m telescope of the Special 
Astrophysical Observatory of the Russian Academy of 
Sciences indicate that the actual period of the outer 
orbit in this triple is twice longer and hence the system 
must be dynamically stable. So far, no systems with 
main-sequence components have been found that could 
be securely classified as dynamically unstable. 

Known candidate objects with low orbit hierarchy, 
which therefore should be viewed as possible dynami- 
cally unstable multiple systems, include the nearby (d 
«16 pc) triple star GJ 795 (HD 196795 = Hip 101955, a 
= 20' l 39 m 38 s , 6 = +04°58'19", epoch 2000.0). Its inte- 
grated spectral type corresponds to that of a K5V star. For 
decades, G J 795 has been known as the visual pair Kui 99 
dKuiper, 193"i| with a period of 40 years. During his spec- 
troscopic survey of visual binaries with the CORAVEL 
radial- velocity scanner, Duquennoy (Duqucnnoy, 1987) 
found GJ 795 to contain a hitherto unknown spectroscopic 
subsystem with very low amplitude of radial- velocity vari- 
ations. He concluded that the companion is bound to the 
main component of the binary and computed a prelimi- 
nary model of the system by combining photometric and 
spectroscopic data with the computed orbital elements. 
At the same time, the inclination of the inner orbit of 
the triple was estimated based on the assumed component 
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Table 1. Log of speckle observations 



GJ795 



Date 


B 

H 


N 


A/AA 




arcsec 




nm 


1998.7741 


1 


700 


610/20 


1999.8206 


2 


1500 


610/20 


2000.8752 


1.5 


1000 


600/30 


2001.7522 


2 


1500 


545/30 




2 


1500 


850/75 


2002.7986 


3-5 


1500 


600/30 


2003.9272 


1 


2000 


545/30 




1 


2000 


700/30 




1 


2000 


800/110 


2004.8232 


1.5 


2000 


600/30 



masses exclusively. The strongly inclined outer orbit also 
remained highly uncertain. To refine the pattern of com- 
ponent motions in GJ 795, this system was put in 1998 into 
the list of program stars for speckle interferometric obser- 
vations with the 6 m telescope of the Special Astrophysical 
Observatory of the Russian Academy of Sciences. 

In this paper, we report the results of our speckle in- 
terferometric observations of the relative measurements of 
components of GJ 795 and their differential photometry, 
and determine the parameters of orbital motion of stars 
and their dynamic masses. In conclusion, we discuss the 
dynamical stability of the system. 

2. OBSERVATIONS AND DATA ANALYSIS 

We performed speckle interferometric observations of 
GJ 795 with the new facility mounted on the 6 m telescope 
of the Special Astrophysical Observatory of the Russian 
Academy of Sciences (Maximo v et al., 2003[ ). Its detector 
consists of a fast 1280 x 1024 Sony ICX085 CCD combined 
with a three-camera image-tube converter with electro- 
static focusing. We recorded speckle interferograms in the 
visible part of the spectrum with exposures ranging from 
5 to 20 milliseconds. Table [1] lists the log of observations, 
which gives for each measurement the date of observation 
(as a fraction of Besselian year); seeing (3 in arcsec; the 
number of speckle interferograms in the series; filter pa- 
rameter A/AA in nm, where A and AA are the central 
wavelength and half-bandwidth, respectively. We deter- 
mined the relative component positions and magnitude 
differences from the power spectra of the speckle inter- 
ferograms averaged over the series jBalega et al., 20 02a 
The accuracy of the measured position parameters is 0.2- 
4.0° and 1-4 milliarcsec for position angle and angular 
separation, respectively. Measurement errors depend on a 
number of parameters: component separation, magnitude 
difference, and seeing [3. The accuracy of the magnitude 
differences inferred from reconstructed power spectra de- 
pends on the same parameters. For objects in the my=8- 
10 magnitude interval it usually varies from 0.05 to 0.2. 
We use bispectral analysis of the interferogram series to 



Fig. 1. The 610/20 nm image of GJ795 (1998.77) recon- 
structed using bispectral analysis. Artefacts surrounding 
the point sources are due to various types. North is at the 
top and East on the left. 



perform complete reconstruction — including that of mod- 
ulus and phase — of the images ( jLohmann et al.,~1 983; 



Weigelt, 1977). Figure [T] shows the reconstructed image of 
the triple star GJ 795 based on observations made in 1998. 

3. ABSOLUTE MAGNITUDES AND SPECTRAL 
TYPES 

As we already mentioned in the Introduction, Duquennoy 
( |Duquennoy, 1987| carried out a detailed study of the 
radial-velocity variations of primary component A of the 
visual binary KUI 99 with the CORAVEL radial- velocity 
scanner. He detected no traces of the fainter star B in 
the spectra. On two nights in 1985, significant varia- 
tions were observed in the profile of the correlation min- 
imum that were due to the contribution of component 
Ab to the total flux. The resulting radial-velocity curve 
was used to determine the orbit of inner binary Aab 
with a period of P=920.2 days and eccentricity e=0.747. 
The preliminary model of the system, which included all 
three components, assumed a total mass and parallax 
of EXab=1-62±0.27A4 q and 7r=64±5 mas, respectively. 
This model also made use of the empirical "mass - lu- 
minosity" relation for K6V (Aab) and K9V(B) type stars 
and a highly uncertain visual orbit of the outer pair AB 
|Baize, 1981] 

According to Hipparcos QPerryman, 1997| ) data, the 
trigonometric parallax of GJ 795 differs significantly 
from the above value ( 717^=53. 82±2. 21 mas). However, 
Hipparcos trigonometric parallaxes for binary and mul- 
tiple stars are known to be potentially fraught with ex- 
tra errors due to wrong correction of component or- 
bital motions in the process of Hipparcos data reduc- 
tion flBalega et al., 2002b[ |Shatskii fc Tokovinin, 19981 ). 
Soderhelm (jSodcrhjclm, 1999) corrected the parallax 
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Table 2. Differential speckle interferometry of GJ 795 



Date 


Comp. 


Am 




A/AA 


R,6f6r6nc6 


BY 


vector 






nm 




1998 7741 


Aa-Ab 


1.09 


0.05 


610/20 


RfllpfTA pi" Fll 




Aa-B 


0.88 


0.05 




2002a 


1999 8206 


Aa-Ab 


1.14 


0.03 


610/90 


J-JdrlCtiCL CLi Oil. 




Aa-B 


0.94 


0.03 




2004 


2000.8752 


Aa-Ab 


1.30 


0.06 


600/30 


Balega et al. 




Aa-B 


1.02 


0.06 




2006 


2001.7522 


Aa-Ab 


1.35 


0.06 


545/30 


Balega et al. 




Aa-B 


1.11 


0.06 




2006 




Aa-Ab 


0.92 


0.06 


850/75 


Balega et al. 




Aa-B 


0.68 


0.06 




2006 


2004.8232 


Aa-Ab 


1.42 


0.05 


600/30 


This paper 




Aa-B 


1.27 


0.05 







for the effect of the orbital motion of the pair AB: 
7r^ ip =58.8±2.1 mas. The corrected Hipparcos parallax 
agrees within the quoted errors with that given by 
Duquennoy QDuquennoy, 1987] ). 

We performed differential speckle photometry of the 
system with the 6 m telescope of the Special Astrophysical 
Observatory of the Russian Academy of Sciences and list 
the results in Table [2 We set the V-band magnitude dif- 
ferences equal to AmAaAb = 1-35 ± 0.06 and ArriAaB = 
1.11±0.06, respectively, implying, given the corrected par- 
allax, the absolute mag nitudes of My a = 7.31 ± 0.08, 
M$ b = 8.66 ± 0.10, and My = 8.42 ± 0.10, respec- 
tively. These absolute magnitudes correspond to the spec- 
tral types of SpAa ~K5, SpAb ~K9, and Sps ~K8, re- 
spectively. The above spectral types agree well with the 
color index B — V=1.24 ( |Strassmeier et al, 2000] ). The 
space velocity components (U, V, W)=(-75. 5,-19. 7,-42. 3) 
( Strassm eier et al., 2000 1 and low emission level in the H 
and K Ca II lines ( Gray et al., 2003] ) imply that this star 
should be classified as a Galactic-disk object with the age 
of 2-3 Gyr. 



4. ORBITAL PARAMETERS AND TOTAL 
MASSES 

The motion of components in a triple system can be sub- 
divided into two components: the motion about the center 
of mass of the inner binary and the motion of the outer 
component about the common center of mass. The orbit 
of the outer system AB was computed by a number of au- 
thors from visual micrometric measurements ( |Baize, 1 981 
|Heintz, 19 84). Its main parameters are: a period P of 
about 40 years; small eccentricity e, and high inclina- 
tion with respect to the sky plane (i «85°). Soderhjelm 
(Soderhjcl m, 1999| refined the orbital elements of the pair 
AB by combining the data of ground-based observations 
with Hipparcos astrometry. He inferred a total mass of 
£.Mab=2.26±0.36.M©, which exceeds significantly the 
mass estimated by Duquennoy ( |Duquennoy, 1987] ) based 
on spectroscopic and visual data. Despite the use of 



Hipparcos astrometry, the outer orbit remains uncertain, 
mostly because of its high inclination. 

We determined the preliminary orbital parameters of 
the subsystems of GJ 795 with allowance for new speckle 
intcrferomeric measurements based on the Fourier trans- 
form of equations of motion (Monet, 1979). At the next 
stage, we refined the orbital elements via differential cor- 
rection based on the least squares method (see comments 
in |Forveille et al., 1999Q . We adopted early visual and sev- 
eral interferometric measurements from the Washington 
Catalog of Binary Stars (http://ad.usno.navy.mil/wds). 
We set the weights of speckle observations made with the 
6 m telescope of the Special Astrophysical Observatory of 
the Russian Academy of Sciences to be 10 times higher 
than those of the visual and interferometric observations 
made with other telescopes. The main reason why we at- 
tributed such low weight factors to earlier data is that 
they did not allow for the binary nature of component A, 
resulting in significant systematic errors. 

We computed the orbit of the inner binary Aab based 
exclusively on the data of interferometric observations 
made with the 6 m telescope of the Special Astrophysical 
Observatory of the Russian Academy of Sciences. Seven 
measurements span over about 2.5 periods along the ap- 
parent ellipse of the orbit (Fig. [2]). An attempt to use ra- 
dial velocities from |Duquennoy, 1 987 to construct a com- 
bined orbit resulted in increased errors of the inferred or- 
bital elements. This is due to the combined effect of small 
radial velocities of Aab system, small inclination of the in- 
ner orbit to the sky plane, and systematic errors in radial- 
velocity measurements due to the influence of the distant 
component B. 

To convert the motion of component B to the center- 
of-mass frame of the inner binary Aab, we set the mass 
ratio equal to qi n ~0.8 based on the empirical "mass - 
My" relation of |Henry k. McCarthy, 1993| and absolute 
magnitudes M^=7.3 and M^ b =8.7. We give the el- 
ements of the outer and inner orbits in Table [31 We 
list all position measurements made with the 6 m tele- 
scope of the Special Astrophysical Observatory of the 
Russian Academy of Sciences and the corresponding resid- 
uals in Table [U The orbital parameters of the inner bi- 
nary Aab agree well with those of the spectroscopic or- 
bit by Duquennoy ( |Duquennoy, 1987] ) and those of the 
outer binary AB, with the refined orbit of Sodcrhclm 
(Soderhjclm, 1999). It is remarkable that the preliminary 
estimates of the inclination and semimajor axis of the 
inner orbit obtained by Duquennoy (Duquennoy, 1987) 
based on published empirical "mass - luminosity" rela- 
tions coincided exactly with their true values inferred from 
the results of our interferometry. 

Let us now determine the angle <fi between the orbital 
planes: 

cos0 = cosi out cosi m + s\iii out sinz m cos(tt out - fi m ), (1) 

where i out and it n are the tilt angles of the orbit of the 
outer and inner binaries with respect to the sky plane, re- 
spectively; Vtout and rii n are the longitudes of the ascend- 
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Fig. 2. Relative ellipses of the orbits of the triple system GJ 795: (a) orbit of binary Aab, (b) orbit of binary AB. The 
filled circles show the speckle interferometric observations made with the 6 m telescope of the Special Astrophysical 
Observatory of the Russian Academy of Sciences; the solid line shows the position of the periastron, and the dotted- 
and-dashcd line, the line of nodes. The radius of the dashed circle is equal to 20 mas. Position parameters of component 
B are converted to the center-of-mass frame of binary Aab. 



Table 3. Parameters of the inner and outer orbits in the 
triple system G J 795 





Aab 


AB 


P, years 


2.51±0.01 


39.4±0.2 


T 


2000.55±0.01 


1975.0±0.3 


e 


0.620±0.006 


0.06±0.01 


a, mas 


120±2 


820±30 


i° 


18±3 


86.9±0.1 


n° 


174±11 


128.5±0.2 


o 


87±11 


212±2 


ere 


0.5 


0.8 


o p 


1 


1 



ing node of the outer and inner binaries, respectively. We 
now use the angle values from Table [3] to obtain </>=74°. 

The dynamic mass of the inner binaries as inferred 
from the orbital parameters and the corrected Hipparcos 
parallax of 7r^ p =58.8±2.1 mas ( |Soderhjelm, 1999D is 
equal to SA^A a 6=1.28±O.15j\40. The mass of the entire 
system GJ 795 computed using the orbital parameters of 
binary AB is equal to J2 .Mab=1.69±0.27./W©. According 
to Lang ( |Lang, 19 92), the individual masses of stars in 
the system as inferred from their absolute magnitudes are 
equal to: A^ Aa =O.67X , M A b=0-57M Q , 7W B =O.54X . 
The above estimates imply a mass of EA^Aafc=l-24A40 
for Aab binary and a total mass of ^ Mab — ^-T^AAq for 
the entire system, which are consistent within the quoted 
errors with the total masses inferred using orbital param- 
eters. 



5. HIERARCHY OF ORBITS AND THE KOZAI 
MECHANISM OF OSCILLATIONS 

The ratio of the orbital periods of the components of the 
system is equal to P ou t/ fm=15.7, and hence the system 
is weakly hierarchical and its components form a single 
gravitationally bound system. Let us now try to estimate 
its dynamical stability using empirical criteria and crite- 
ria based on numerical simulations. According to one of 
such criteria suggested by Tokovinin ( |Tokovinin, 2004 1, 
the system is stable if the following inequality is satisfied: 



T 



Print ( 1 



p. 

- 1 1.1 



(2) 



where P ou t and Pi n are the orbital periods of the outer and 
inner binaries, respectively; e ut, the eccentricity of the 
orbit of the outer binary, and T c is the critical instability 
value, which is equal to 5. The stability parameter for 
GJ 795 is equal to T«13 and hence the system is stable. 

Based on numerical simulations, Harrington 
(Harring ton, 1977| proposed the following stability 
parameter for triple systems, which depends on the ratios 
of the semimajor axes of the outer and inner orbits, a ou t 
and a,i n , and the eccentricity of the outer orbit: 



r = > r c . 



(3) 



Like in the case of empirical estimate, the value of this 
parameter, F=6.42, exceeds the critical level of F c =5.46. 
However, these criteria should be used with much caution 
when applied to systems with orthogonal orbits. 

The semimajor axes of the outer and inner binaries 
are equal to a out «14 AU and <ij n «2 AU, respectively. 
The apoastron distance of the inner binary is equal to 
«3.2 AU. The eccentricity of the outer orbit is close to 
zero and therefore components of the triple are never at 
comparable distances from each other. 
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Theoretical studies on the dynamics of multiple sys- 
tems show that in the case of large angles between the 
orbital planes the inner and outer binaries ( |Kozai, 1962[ ) 
exchange angular momentum. This mechanism triggers 
periodic variations (Kozai oscillations) of the eccentric- 
ity of the inner orbit, e; n , and angle <f> between the orbital 
planes. The quantity (1 — ef n ) cos 2 4> = const remains con- 
stant in this process. The following formula gives the pe- 
riod of Kozai cycle: 

Pkozai ~ Pout/Ptnil - e out f 2 . (4) 

The period of Kozai oscillations for the triple star 
GJ 795 is equal to only 560 years. We may try to directly 
observe Kozai oscillations of the orbital parameters of bi- 
nary Aab over several years of interferometric observa- 
tions. 

6. CONSLUSIONS 

We use speckle interferometric observations made in 1998- 
2004 with the 6 m telescope of the Special Astrophysical 
Observatory of the Russian Academy of Sciences to com- 
pute accurate visual orbits for the outer and inner binaries 
of the triple star GJ795. This triple system belongs to 
the disk component of the Galaxy and is 2-3 Gyr old. 
Differential photometry of the components of this sys- 
tem made it possible to construct a complete model of 
GJ 795, which agrees well with modern empirical and the- 
oretical relations. The absolute magnitudes of the compo- 
nents are equal to M^ a =7.31±0.08, M^ b =8.66±0.10, and 
My =8.42±0.10 and they correspond to the spectral types 
of SpAa ~K5, SpAb ~K9, and Sps ~K8, respectively. 

The orbital periods are equal to 2.51 and 39.4 years 
for the inner and outer binaries, respectively. The an- 
gle between the planes of the inner and outer or- 
bits is equal to 0=74°. The total dynamical masses 
£.M Aab =1.28±O.15.M and Y.-^ab = 1.69±O.27A4 
are consistent with the estimated spectral types of the 
components. 

We use the available empirical and theoretical stability 
criteria to conclude that GJ 795 is a gravitationally bound 
stable hierarchical system. For objects of this type, the 
Kozai mechanism should be efficient, which causes oscilla- 
tions of the orbital eccentricities and the angle between the 
orbital planes. During its lifetime, the triple star GJ795, 
whose Kozai period is equal to Pkozai ~560 years, must 
have undergone ~ 10 6 such periodic perturbations. 
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Table 4. Position parameters and residuals of the measurements of the triple : 



l GJ 795 
system GJ 795 



Subsystem 


Epoch 


e 


p 


(O - C)o 


(o-c) p 






degrees 


mas 


degrees 


mas 


Aab 


1998.7741 


55.0 


161 


0.3 


-2 




1999.8206 


105.7 


164 


-0.5 


1 




2000.8752 


20.8 


108 


0.4 


1 




2001.7522 


78.2 


185 


-0.2 







2001.7522 


78.5 


186 


0.5 


1 




2002.7986 


150.3 


98 


0.9 


2 




2003.9272 


61.9 


172 


-0.3 


-1 




2003.9272 


61.9 


174 


-0.3 


1 




2003.9272 


61.7 


174 


-0.5 


1 




2004.8232 


105.2 


164 


-0.1 





AB 


1998.7741 


135.4 


358 


0.2 


2 




1999.8206 


139.1 


234 


-0.2 


1 




2000.8752 


153.2 


108 


0.0 


1 




2001.7522 


238.5 


48 


2.3 


1 




2001.7522 


236.8 


45 


0.6 


-2 




2002.7986 


292.5 


153 


0.6 







2003.9272 


300.3 


286 


0.1 


-3 




2003.9272 


300.0 


289 


-0.2 


1 




2003.9272 


299.8 


289 


-0.4 


1 




2004.8232 


303.1 


390 


0.3 


-1 



